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ABSTRACT 

We extend the analysis of the MAXIMA-1 cosmic microwave background (CMB) data to smaller 
angular scales. MAXIMA, a bolometric balloon experiment, mapped a 124 deg 2 region of the sky with 
10' resolution at frequencies of 150, 240 and 410 GHz during its first flight. The original analysis, 
which covered the multipole range 36 < I < 785, is extended to I = 1235 using data from three 150 GHz 
photometers in the fully cross-linked central 60 deg 2 of the map. The main improvement over the original 
analysis is the use of 3' square pixels in the calculation of the map. The new analysis is consistent with 
the original for t < 785. For I > 785, where inflationary models predict a third acoustic peak, the new 
analysis shows power with an amplitude of 56 ± 7 /iK at £ 850 in excess to the average power of 42 ± 3 
//K in the range 441 < i < 785. 

Subject headings: cosmic microwave background - cosmology: observations 



1. INTRODUCTION 

The power spectrum of temperature fluctuations in 
the cosmic microwave background (CMB) has recently 
been measured over a broad range of angular scales (de 
Bernardis et al. 2000; Hanany et al. 2000). The results are 
being used to discriminate between cosmological models 
and to determine cosmological parameters (Lange et al. 
2000; Balbi et al. 2000; Dodelson & Knox 2000; Tegmark 
k Zaldarriaga 2000; Pierpaoli, Scott, & White 2000; Jaffe 
et al. 2000). There is strong evidence for a first peak in the 
power spectrum at I ~ 200 which is consistent with a flat 
universe in adiabatic inflationary models with cold dark 
matter. Fluctuations were also detected in the multipole 
range 400 < I < 785 at smaller scales than those of the 
first peak. 

MAXIMA is a balloon-borne experiment optimized to 
map the CMB anisotropy over a broad range of angular 
scales from 10' to 5°. The first flight, MAXIMA-1, ob- 
served 124 square degrees of the sky. We have published an 



analysis of this data over the spherical harmonic multipole 
range 36 < I < 785 (Hanany et al. 2000). In this paper, 
we extend the original analysis to £ — 1235 using the data 
from three 150 GHz photometers. These measurements 
are a powerful discriminant between cosmological models 
and allow parameters to be determined with better accu- 
racy. A companion paper Stompor et al. (2001) discusses 
the cosmological significance of these results. 

The MAXIMA instrument is described in detail by Lee 
et al. (1999) and Hanany et al. (2000). In sections 2-4, 
we describe the analysis of the telescope beams, the point- 
ing reconstruction, and the estimation of instrument noise, 
which are especially important for the high-£ region of the 
power spectrum. In section 5, we present the map and 
the angular power spectrum. In sections 6 and 7, we show 
that foreground contributions to the data are negligible 
and describe tests for systematic errors. 

2. TELESCOPE BEAMS 
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For a total-power radiometer such as MAXIMA, a finite 
telescope beam size introduces an exponential reduction 
in sensitivity at high-£ Extending the MAXIMA power 
spectrum to higher t values thus depends on an accurate 
knowledge of the beam patterns. Each of the 16 elements 
of the MAXIMA bolometer array were designed to have a 
FWHM Gaussian beam size of 10', which reduces power 
by a factor of five at I ~ 1000. The beams were mea- 
sured in flight using observations of Jupiter, which has an 
angular diameter of 0.5'. Observations were carried out 
at fixed elevation with the primary mirror scanning ~ 4° 
peak-peak in azimuth at 0.45 Hz and the gondola azimuth 
fixed. Sky rotation gave a raster map of Jupiter with more 
than 200 mirror scans across each array element. The re- 
constructed beam profiles have a peak signal-to-noise ratio 
of more than 1000. The beams are close to circular near 
the peak, but have a moderate asymmetry in the wings 
(see Fig. 9 in Wu et al. (2000)). 

Wu ct al. (2000) have developed a formalism for cal- 
culating the effects of asymmetric beams. We find that 
data measured with slightly asymmetric beams can be an- 
alyzed by using an effective symmetric beam. This formal 
treatment has been tested using simulated observations of 
a CMB sky with the observation strategy and measured 
beams for MAXIMA- 1. Because of the high degree of sym- 
metry near the beam centers, the estimated error in Ce 
introduced by this technique is less than 1% at I = 1000. 

The uncertainty in the measured beam profiles is a 
source of error at high I. There is a 5% uncertainty in the 
measurement of the long and short axis FWHM, which is 
largely due to uncertainty in the baseline subtraction for 
each of the one-dimensional Jupiter scans. There are mi- 
nor contributions to this error from uncertainties in point- 
ing reconstruction and the measured bolometer time con- 
stant. This error gives an error in Ce that increases with 
I (Wu et al. 2000). The error is 6% at I = 500 and 17% at 
I = 1000. The errors in the t bins are correlated, so the 
power spectrum will be biased either up or down with an 
amplitude that increases with I. 

3. POINTING RECONSTRUCTION 

The MAXIMA experiment observes each sky pixel re- 
peatedly to achieve an adequate signal-to-noise ratio. Er- 
rors in telescope pointing reconstruction smear the signal 
between pixels. 

There are four temporal modulations of the beam in 
the MAXIMA experiment. The primary mirror scans in 
azimuth by ~ 4° peak-peak at ~ 0.45 Hz, the gondola 
scans in azimuth by ~ 8° peak-peak at ~ 20 mHz, sky ro- 
tation is used to scan in right ascension, and we map the 
same sky region twice during the flight at different eleva- 
tions to achieve cross linking of pixels. Information from 
several sensors is combined to give the telescope pointing 
as a function of time, and each sensor contributes to the 
total pointing error. 

The dominant pointing error comes from the primary 
mirror motion, as the mirror is rotated left or right. The 
elevation angle of the beams increases as the square of the 
azimuth scan angle. This elevation change can approach 
several arc minutes. We measure the centering of the pri- 
mary mirror motion during flight with a ~ 1° uncertainty, 
which gives a 0.8' elevation uncertainty. 



A CCD camera mounted to the inner frame of the gon- 
dola measures the attitude of the telescope by monitoring 
stars. A dedicated processor calculates the position of the 
two brightest stars in the field of the camera. By compar- 
ing the reconstructed and known separations of the stars, 
we estimate a 0.45' rms error in gondola orientation. 

The total pointing reconstruction error is 0.95' rms, in- 
cluding the described contributions from the primary mir- 
ror position and the gondola attitude errors and several ad- 
ditional minor contributions. Although some of the point- 
ing errors are systematic rather than random, the multiple 
modulations in MAXIMA'S observation strategy tend to 
randomize the errors when projected on the sky. There- 
fore, we model the pointing error as a simple Gaussian 
distributed error. 

We simulate an observation with such a Gaussian point- 
ing reconstruction error to assess the effect on the mea- 
sured power spectrum. As expected, the main effect is a 
systematic underestimate of power at high L In principle, 
this bias is predictable and can be accounted for accu- 
rately. However, since the pointing error is not perfectly 
Gaussian, we conservatively assign a symmetric random 
error equal to the size of the bias which is SCe/Ce ~ 10% 
at I = 1000. The Ce error for each I bin due to pointing 
reconstruction error is given in Table 1. As with the beam 
shape uncertainty, this fully correlated error has the ef- 
fect of increasing or decreasing power monotonically with 
increasing I. 

4. NOISE ESTIMATION 

The data contain both the signal from the sky and noise 
from the detector system. The procedures for map mak- 
ing and power spectrum estimation require an accurate 
estimate of the detector system noise to correctly estimate 
the signal. The accuracy of the noise estimation is most 
critical in the high-^ portion of the power spectrum, where 
the data become increasingly noise dominated. This trend 
occurs because Ce falls roughly as £ 2 and because the tele- 
scope beam reduces the signal at high I. 

As described in (Hanany et al. 2000), we observe a spuri- 
ous primary-mirror synchronous signal which complicates 
the noise estimation because it is larger than the rms am- 
plitude of the noise. We use two methods to estimate the 
detector noise. In the first method, we assume that noise 
dominates the signal at all frequencies other than the pri- 
mary mirror modulation frequency and its harmonics. The 
noise at those frequencies is calculated by interpolating 
data from the power spectrum at nearby frequencies. In 
the second method, we use the iterative noise estimation 
procedure of Ferreira & Jaffe (2000) . Here, the noise and 
signal are both estimated iteratively in the map-making 
stage. We iterate the noise using a map with 8' pixels to 
limit computation time. This resolution should be suffi- 
cient, since the signal power decreases at smaller angular 
scales. 

We make 3'-resolution maps using the two noise esti- 
mates. We calculate the power spectrum for <~ 23, 000 
pixels in the "central" region of the maps, and we find 
that the difference between the two spectra is much less 
than the statistical error for all t bins. This central region 
covers roughly half the area of the map (60 deg 2 ) where 
the observations are fully cross linked, the sampling is the 
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most uniform, and the signal-to-noise ratio is highest. In 
this way, we reduce both the time required to search for 
systematic errors and the computation time. This proce- 
dure increases statistical error bars by < 10% for £ > 500 
and by < 3% at £ = 1000. 

The power spectrum of the difference between maps 
from two photometers in the array is a stringent test of 
the noise estimation process. A mis-estimate of the white 
noise level of the detector system gives a power spectrum 
that diverges from zero as £ 2 . Correlation in the noise be- 
tween detectors can give the same behavior. In a separate 
analysis we do not find significant correlation between the 
time streams in the 0.1 - 30 Hz frequency range used in the 
analysis aside from that caused by the spurious chopper- 
synchronous signal. This signal is removed during the map 
making stage resulting in maps without correlations. We 
also find that histograms of the temperatures in the 8' 
difference maps are consistent with the distributions ex- 
pected for no noise correlations at a Kolmogorov-Smirnov 
significance level larger than 10%. As discussed in section 
7, the power spectra obtained using difference maps are 
consistent with zero, with no apparent systematic diver- 
gence at high £. 

5. MAP AND ANGULAR POWER SPECTRUM 

The map-making and power spectrum analysis proce- 
dures are similar to those described by Hanany et al. 
(2000). In that analysis we used 5' x 5' pixels to limit com- 
putation time, and we limited the analysis to an £ range 
where such pixelization does not bias the results. In prin- 
ciple, the smearing effect of larger pixels can be accounted 
for by deconvolving a pixel window function. In practice, 
however, this procedure may introduce a systematic bias 
in the power spectrum due to the uneven spatial distri- 
bution of samples in each pixel. We choose a 3' square 
pixel for the present analysis, since this pixel size does not 
compromise the resolution of the map (Wu et al. 2000). 

In this paper, we report on the analysis of data from 
three 150 GHz photometers using a 3' pixelized map. 
Hanany et al. (2000) analyzed the data from three 150 
GHz photometers and a 240 GHz photometer using 5' pix- 
elization. The data from all four photometers passed all 
tests for consistency up to £ = 785. In this work with 3' 
pixels, we exclude the 240 GHz data because it does not 
pass consistency tests above £ — 785. The power spec- 
trum of difference maps between the 240 GHz data and 
any 150 GHz data are inconsistent with zero in this high 
£ region. One possible reason for this discrepancy is that 
the spurious primary-mirror synchronous signal is larger 
by a factor > 2 in the 240 GHz data compared to the 150 
GHz data. 

The raw data preparation is described in Hanany et al. 
(2000). We remove data during calibration events, and we 
also remove short sections of data which deviate by > 4 a 
such as cosmic ray events and telemetry drop-outs. We 
deconvolve the transfer functions of the bolometers and 
readout electronics and estimate the noise power spectrum 
from sections of the time stream that had no long gaps 
(Stompor et al., in preparation). We marginalize over fre- 
quencies lower than 0.1 Hz and higher than 30 Hz, where 
we do not expect appreciable optical signals. 

The data are calibrated using observations of the CMB 



dipole, and periodic optical stimulator events are used to 
account for a small drift in calibration during the flight 
(Hanany et al. 2000). The calibrated time stream data, 
the pointing solution, and the estimate of the noise power 
spectrum are combined to produce a maximum likelihood 
pixelized map of temperature anisotropy and a pixel-pixel 
noise correlation matrix for each photometer (Wright 1996; 
Tegmark 1997; Bond et al. 1999) using the MADCAP soft- 
ware package (Borrill 1999). 

A combined maximum-likelihood temperature 
anisotropy map is produced by adding individual maps 
with a weight inversely proportional to their noise corre- 
lation matrices. This map, shown in Figure 1, contains 
~ 40, 000 3' square pixels. We assign a calibration uncer- 
tainty of 4% to the magnitude of temperature fluctuations 
in the combined map (Hanany et al. 2000). 

We calculate the angular power spectrum Ce of the cen- 
tral region of the combined map using the MADCAP (Bor- 
rill 1999) implementation of the maximization of the like- 
lihood function following Bond, Jaffe & Knox (1998). The 
MADCAP implementation assumes that the beam shape 
has axial symmetry. We produce an effective beam for the 
analysis of the combined map by noise- weighted averaging 
the individual beams. As discussed in section 2, we fol- 
low the procedure of Wu et al. (2000) to find a symmetric 
approximation for the effective beam. We apply no cor- 
rections for the negligible (< 1% at £ = 1000) smoothing 
of the 3' pixels. 

We calculate the power spectrum of the temperature 
fluctuations using 15 bins in £ over the range 3 < £ < 1800 
assuming a constant £(£ + l)Cg band power in each bin, 
and marginalizing over the CMB monopole and dipole. 

The benefits of only using the central region of the map 
for power spectrum analysis are discussed in section 4. 
The penalty is an increased error that falls with increas- 
ing £ due to the rise in sample variance. To minimize this 
effect, we create a composite power spectrum which uses 
the £ < 335 points from the original Hanany et al. (2000) 
analysis. We choose this transition point to limit the max- 
imum error bar increase to ~ 20%. For bins with £ < 335, 
we marginalize over all the higher £ bins from Hanany et 
al. (2000). For the bins with £ > 335 obtained using the 
central region, we marginalize over the bins at £ < 335 
and £ > 1235. For both £ regions, we diagonalize the £-bin 
correlation matrix using a variant of techniques discussed 
in Bond, Jaffe & Knox (1998). The correlations between 
the dominant bin and adjacent bins are typically less than 
10%. 

Table 1 lists the dominant bins, the Ci estimates, and 
the AT = \J £{£ + l)Ce/2TT estimates for the correspond- 
ing uncorrelated linear combinations of bins. We quote 
la errors on the d estimates assuming 68% confidence 
intervals using the offset log-normal distribution model of 
Bond, Jaffe & Knox (2000). These statistical errors do 
not include three additional independent sources of sys- 
tematic uncertainty, which arc fully correlated between 
£ bins. The la calibration error is a constant 8% of 
£{£ + 1)Ci/2tt. The ^-dependent errors due to the beam- 
diameter uncertainty and pointing reconstruction uncer- 
tainty are given in Table 1. Information on the shape 
of the bin-power likelihood functions and window func- 
tions will be made available on the MAXIMA web site 
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Table 1 

Uncorrelated Power Spectrum from the MAXIMA- 1 map. Data for £ < 335 are from the power spectrum of the 
full map with 5' pixels, and data for £ > 335 are from the power spectrum of the central region of the map 

with 3' pixels. Errors are the 68% integrated pro bability of t he offset log-normal likelihood functions with 
a constant prior in either £(£ + l)Ci/2% or AT = \J l{£ + 1)Ci/2tt for their respective columns. Here [£min,£max] 

GIVES THE RANGES OF THE DOMINANT BINS. THE CORRELATED BEAM ERRORS ARE DESCRIBED IN THE TEXT. 



(http: / / cfpa.bcrkeley.edu/maxima) . 

The top panel of Figure 2 shows the maximum likeli- 
hood power spectrum and an inflationary adiabatic model 
that best fits the Hanany et al. (2000) MAXIMA-1 and 
COBE/DMR power spectra. The results for £ < 785 are 
consistent with those in Hanany et al. (2000). The model 
has (n b , n cdm , A , n, h) = (0.07, 0.61, 0.23, 1, 0.60) (Balbi 
et al. 2000). The \ 2 for this model is 37 for all 41 data 
points. If we only use the 13 data points of MAXIMA-1 
we obtain x 2 = 5 f° r tms model. In both cases, we are 
fitting with an 7 parameter inflationary model. 

6. FOREGROUNDS 

Foreground sources include: emission from the earth, 
the atmosphere, and galactic dust, free-free and syn- 
chrotron radiation, point sources, and scattering due to 
the Sunyaev-Zel'dovich effect. The analysis of Hanany ct 
al. (2000) for £ < 785 is applicable for sources that are 
expected to be less strong at small angular scales. In that 
work, we argue that the detected signal is inconsistent with 
an atmospheric or ground-based origin due to its temporal 
stability. We can detect dust contamination in our maps, 
but at a negligible level of 2.3 /xK at 150 GHz. Estimates 
of bremsstrahlung and synchrotron radiation (Bouchet & 
Gispert 1999) predict contributions of less than 1 /iK at 
150 and 240 GHz. The ratio of detected power between 
150 and 240 GHz is consistent with the CMB and incon- 
sistent with dust, synchrotron, and free- free emission. 

Point sources are the foreground of concern for the high- 
£ region data presented in this paper. A catalog search 
(Sokasian, Gawiser, & Smoot 2000; Gawiser & Smoot 
1997) yielded no detectable radio or infra-red sources in 
any of the frequency bands. The "pessimistic" model of 
Tegmark et al. (2000) predicts ~ 20 fiK rms point source 
contamination at £ = 1000, which is small compared to 
our statistical errors. The shape of the power spectrum is 



inconsistent with a point source origin, which would give a 
rising spectrum with increasing £. Instead, the power spec- 
trum decreases above £ — 850 and is consistent with zero. 
A broad bin covering 1235 < £ < 1800 (not tabulated) is 
consistent with zero at the 1 a level. 

7. TESTS FOR SYSTEMATIC ERRORS 

In the original data analysis of Hanany et al. (2000), we 
found no evidence for uncorrected systematic errors over 
the range 35 < £ < 785 after performing a suite of tests. 
Although maps with 8' and 10' resolution give power spec- 
tra with a systematic bias at high £, they are still useful 
for systematic error tests especially in the case where the 
result is a null power spectrum. We confirm that the fol- 
lowing combinations of data produce a spectrum with no 
signal over the range 35 < £ < 1235: (1) a dark bolome- 
ter, (2) the data from the 410 GHz photometer, (3) the 
difference between the overlapping part of the combined 
map from the CMB-1 and -2 scans, and (4) the differences 
between the maps produced by different photometers. We 
repeated the fourth test on the central region of the map 
using 3' pixels. We calculated the difference of one and 
the sum of the other two as shown in Fig. 2 

8. DISCUSSION 

We have presented a measurement of the angular power 
spectrum of the CMB over a range of angular scales cor- 
responding to the multipole range 36 < £ < 1235, which 
is the largest yet reported by a single experiment. The 
data for £ > 400 can be fit to a flat line of amplitude 1860 
fiK 2 with a x 2 ~ 5 for 7 d.o.f. However, if we focus on 
the bin at £ ~ 850, we find that its amplitude of 56 ± 7 
[iK is in excess to the average power of 42 ± 3 fiK in the 
range 441 < £ < 785 at a confidence level of ~ 95%. This 
confidence level calculation is presented in the companion 
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paper by Stompor et al. (2001). This excess power is con- 
sistent with the third acoustic peak predicted in this region 
by inflationary adiabatic models. The power spectrum is 
well fit by such a model over the entire range of I. 
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Fig. 1. — A maximum-likelihood map of the CMB anisotropy from MAXIMA-1. The resolution is determined by the 10' FWHM Gaussian 
beam of the telescope. The map is made using data from three 150 GHz photometers and contains ~ 40, 000 3'x3' pixels. We use the central 
region of the map (outlined) to compute the power spectrum presented in Fig. 2. This region is fully cross linked, has the most uniform 
sampling, and has the highest signal-to-noisc ratio. The area of the central region is 60 deg 2 and contains ~ 23, 000 pixels. 
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Fig. 2. — Top panel: Composite angular power spectrum of the CMB anisotropy from the MAXIMA-1 map (filled circles). The points for 
t < 335 are from the power spectrum of the full 5' pixelized map from Hanany et al. (2000), and the points for I > 335 are from the power 
spectrum of the central region of the 3' pixelized map shown in Fig. 1. The error bars are 68% confidence intervals calculated using the offset 
log-normal likelihood functions of Bond, Jaffe & Knox (2000). The solid curve is the best fit (ACDM) inflationary adiabatic cosmology to the 
MAXIMA-1 and COBE/DMR data. The model has (Q b , Cl cdm , Q A , n, h) = (0.1, 0.6, 0.3, 1.08, 0.53) (Balbi et al. 2000). The crosses are the 
estimated power spectrum of the difference between two independent maps, the first given by one of the three photometers and the other from 
the sum of the other two. Bottom panel: A comparison between the Hanany et al. (2000) power spectrum and the power spectrum of the 
central region of the 3' pixelized map shown in Fig. 1. For the 3' central region data, the power spectrum for I < 186 is not well constrained. 



